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Abstract 

We derive model-independent lower bounds on the stress tensor central charge Ct in 
terms of the operator content of a 4-dimensional Conformal Field Theory. More precisely, 
Ct is bounded from below by a universal function of the dimensions of the lowest and second- 
lowest scalars present in the CFT. The method uses the crossing symmetry constraint of the 
4-point function, analyzed by means of the conformal block decomposition. 
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1 Introduction 



Conformal Field Theory (CFT) was born to describe fixed points of renormalization group flows. 
This still remains its main vocation, although it has many other applications as well. In 2D, the 
constraining power of conformal symmetry is tremendous and often leads to an exact solution of 
the theory. In this paper we are concerned with the much less constrained 4D case. Presumably, 
there are lots of interacting 4D CFTs out there, but we don't know much about them. For instance, 
"conformal windows" of gauge theories should provide lots of examples, but even the spectrum of 
these theories remains unknown (except for the chiral ring in the supersymmetric case). 

In absence of an exact solution, it is natural to look for universal constraints, satisfied all over 
the "landscape" of CFTs. Two years ago [T], [2] we found one such constraint, related to a gap in 
the spectrum of operator dimensions. Namely, we examined the maximal possible dimension of 
the lowest-dimension operator appearing in the Operator Product Expansion (OPE) of two scalar 
operators. We found that if one fixes the dimension of external scalars, the lowest-dimension 
operator that appears cannot have a dimension above a certain model-independent bound. 

Last year a different constraint was presented in [3J: it was found that the OPE coefficient of 
three scalars cannot exceed a certain universal 0(1) bound, which depends only on their dimen- 
sions. One can call this bound a universal limit on the interaction strength. 

The above results were obtained by using consistency between OPE and crossing symmetry 
(also known as OPE associativity) as a constraining principle. The principle itself was first pro- 
posed more than 35 years ago by Polyakov [3] , but until our work no general results were obtained 
from it. 

We would like to mention a parallel line of development in 2D CFT, where Hellerman [5] and 
others [6], [7] have also studied constraints on the gap in the spectrum, in terms of the central 
charge. Their main constraining principle is modular invariance, which is limited to 2D, but 
morally not so different from OPE associativity (both are related to the change of foliation when 
quantizing the theory). Also interesting is the role played by all these universal constraints in an 
ambitious program of exploring the space of CFTs initiated by Douglas [8]. 

Coming back to 4D, in this note we will explore the following question: What can we say 
about the central charge of the theory, if we know something about the spectrum of its operator 
dimensions? More precisely, we will assume that the theory contains a scalar operator of a given 
dimension. Under this assumption, we will show that the central charge must be bigger that a 
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certain universal lower bound. This is natural, since central charge "measures" the number of 
degrees of freedom, and by assumption we know that our theory is not trivial. 

In a certain range of scalar dimensions we will be able to show that the central charge is 
necessarily bigger than that of the free scalar ("interacting theory has more degrees of freedom"). 

As we will explain below, the problem of bounding the central charge from below is equivalent 
to the problem of bounding from above the OPE coefficient of the stress tensor in the scalar times 
scalar OPE. The similarity with the problem analyzed in [3] is then clear, and we will be able to 
use the method of that paper. To make connection with the results of [lj,|2j, we will also study 
how the central charge bound improves as a function of the assumed gap in the scalar sector of 
the OPE. 

Everywhere we assume that we are dealing with a unitary theory. In a non-unitary theory, 
central charge may well be zero (or negative) without CFT being trivial, so our question would 
not even make sense. 



2 Formulation of the problem 



We begin by stating precisely our assumptions and goals. 

In an arbitrary unitary CFT in D = 4 spacetime dimensions, we consider the central charge 
Ct, defined as the coefficient in the 2-point function of the stress tensor operator T^: 
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[Sn = 2ix D l 2 /T(D /2)). It is assumed that the stress tensor is normalized canonically, that is 
consistently with the Ward identities (written here schematically for scalars): 



df, {T tlv {x)4>{x 1 ) . . . 4>(x n )) = -y 8(x - Xj) (0(xi) . . . d u (j){xi) . . . 4>(x n )) . 



(2) 



The central charge Ct is an interesting quantity because it provides a certain measure of the 
number of degrees of freedom in the theory. For example, for a free conformal theory of 
scalars, N$ Dirac fermions, and Na vectors, we have [9] 

4 



C T = -iV> + 8iV> + 16N A . 



(3) 
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Moreover, by unitarity Ct > 0, and Ct — corresponds to a trivial theory. It is well known that 
Ct is an imperfect measure since, unlike in 2D, it does not in general decrease along the RG flow 
|10j . The other central charge a, defined in terms of the 4D trace anomaly, fares better in this 
respect [llj . while still remaining imperfect [12]. In this paper, we will stick to Ct since it is the 
one which we are able to constrain. 

Now assume that our theory contains a primary Hermitear^j] scalar operator of a given 
dimension d. Our main goal will be to show that under this assumption, the central charge of the 
theory cannot become arbitrarily small. In other words, there exists a certain universal bound 



where f(d) depends on d but is otherwise model- independent. In this paper we will derive such a 
bound in the interval 1 < d < 2. 

3 Solution strategy 
3.1 Conformal blocks 

We will approach this problem by imposing the constraint of OPE associativity in the 4-point 
function of the operator 0. Consider all primary Hermitean operators appearing in the OPE 

0x0: 



Here in the first line we included the unit operator and all scalar primaries, starting from a certain 
dimension A > 1 and higher. In the second line we have the stress tensor (spin 2, dimension 
4 primary) and possibly higher dimension spin 2 fields. The third line contains all higher spin 
primaries (A > Z + 2 by the unitarity bounds [13j ) . Note that by permutation symmetry of the 
00 state only even spins can appear in this OPE. 

Now, it has been shown by Dolan and Osborn jH] that every primary spin I, dimension A 
operator 0^,; appearing in the 0x0 OPE with a coefficient ca,z gives a contribution to the 4-point 

1 If (ft is not Hermitean, we can consider its real and imaginary parts. 



C T > f(d) > 0, 



(4) 



0x0d1 + Sa + ... (spin 0) 
+ T IMI + ... (spin 2) 
+ higher spins . 



(5) 
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function of of the following form: 
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The functions of the cross-ratios g&.,i(u, v) are called conformal blocks. They can be thought of as 
summing up the contributions of the primary A j and all its descendants to the 4-point function 
of <f), when applying the OPE in the (12) (34) channel. It is nontrivial that such summation can 
be performed in closed form. The clear advantage of the representation (|6| is that it can be used 
at finite point separation, unlike the OPE which is only useful in the coincidence limit Xi — > Xj. 



3.2 Normalizations 



Eq. ^ assumes that both <p and O are unit normalized: 
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and the coefficient ca,i is extracted from the 3-point function 
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At the same time for the stress tensor normalized canonically ([T]),([2]), the 3-point function coeffi 
cient is fixed by the Ward identity 0: 
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These relations determine the coefficient 04,2 appearing in (|6j) and (10) in terms of the central 
charge Ct and the dimension of <p [14J: 

Dd 1 
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Via (|6]), this crucial relation implies that for large Ct, the contribution of the stress tensor to 
the 4-point function of <f> decreases as 1/CV- This is in accord with what happens for example in 
AdS/CFT [15], where Ct ~ N 2 while the stress tensor contribution corresponds to the graviton 
exchange in the bulk, which is 1/N 2 suppressed. Theories without stress tensor (i.e. in which 
c 4,2 = 0) should be viewed as theories with an infinite central charge. One example of such a theory 
is the Gaussian scalar field of dimension d > 1, see [16] for its conformal block decomposition. 

However in this paper we are interested in constraining the opposite limit of small Ct, in which 
the stress tensor contribution increases. As we will see below, such an increase eventually becomes 
inconsistent with crossing symmetry, and this will give us the bound Q. 



3.3 Analytic structure of conformal blocks and crossing symmetry 

Let us now discuss the analytic structure of the conformal blocks. The variable z appearing in 
([7]) may look ad hoc, but in fact it is the 4D analogue of the usual complex variable of the 2D 
CFT, see Fig. [TJ In the Euclidean signature z is complex and z = z*. In this case the conformal 
blocks are real functions, smooth everywhere away from z = and from the (l,+oo) cut along 
the real axis. Since the imaginary part of the hypergeometrics is discontinuous across the cut, 
the conformal blocks have a l/Imz singularity there. Everywhere else on the real axis, and away 
from it, they are regular. 




Figure 1: Using conformal freedom, any configuration of 4 points can be mapped into the one 
shown in this figure, in which 3 points are fixed and one (x^) is moving in a two-plane passing 
through X\ and x 3 . The complex coordinate of x 2 in this plane is precisely the z of (T?p ; while 
z = z* . The conformal blocks are smooth everywhere in the plane except for z = and the shown 
(1, +00) cut along the real axis. 

The asymptotic behavior of conformal blocks as z — > is fixed by the OPE. The singular 
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behavior in the z — > 1 limit, which corresponds to the crossed channel x 2 — > £3, has no simple 
physical meaning. However, the sum over all blocks must be crossing symmetric and thus consis- 
tent with the OPE in the crossed channel as well. Because of the unphysicall singularities, it is 
not immediately clear how to impose the OPE consistency in the crossed channel. In fact for this 
reason Polyakov [I] suggested to use a different type of expansion, into objects he called unitary 
blocks. However at the time the explicit and simple expressions ^ were of course not yet known. 
Armed with these expressions, a different strategy becomes possible. 

Namely, we will study the crossing symmetry condition at finite point separation, which can 
be written as 

G(u,v) _ G{v,u) ^ 



where we used the fact that crossing x\ -H- x% corresponds to the interchange of u and v. Here 
G(u, v) is the sum over all contributing conformal blocks: 



G(u, v) = 1 + ^ c a,i 9a,i(u, ■ 



(15) 



where c\ x are the squares of the OPE coefficients, and we separated the contribution of the unit 
operator. It will be important that in a unitary theory all c^,i are real, so that their squares are 
positive pp. 



r x 4 




Figure 2: This configuration, with 4 points at the vertices of a square, is conformally equivalent 
to the one in Fig. [7] with z = 1/2. 



Instead of going straight to the crossed OPE limit x 2 — > x 3 , we will study (14) around the 
democratic configuration^ when x 2 is at equal distances from x\ and 23. This corresponds to 
z = 1/2. In fact the same configuration can be mapped conformally to 4 operators inserted at 



the vertices of a square (Fig. \2y. Both sides of (14) are regular around z = 1/2. Expanding the 



2 Similarly, Hellerman [5] in his analysis of the modular invariance constraint chose to work at the selfdual inverse 
temperature (3 = 2tt. 
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crossing condition into a two-dimensional power series around this point, we get a infinite number 
of linear equations, which have to be satisfied for some positive coefficients c 2 Al . Which A,/ will 
enter the expansion with nonzero coefficients depends on the CFT. The problem of unphysical 



singularities, brought up by Polyakov jl], is resolved as follows. The LHS of (14) is smooth away 
from the cut along (1, +00), while the RHS away from (—00, 0), since the crossing maps z — > 1 — z. 
Assuming that both sides can be analytically continued from their region of convergence, the cuts 
must cancel when summing over all A and I. In other words, imposing that there is no cut should 
give no additional constraints compared to the ones that we are using, although it may be a 
different and perhaps a useful way to package the same information. 



3.4 Method of linear functionals 



For further discussion let us rewrite Eq. (14) in the equivalent 'sum rule' form: 



1 = X! C A>' F <W(M> ( 16 ) 

j? ( \ v d g A j(u,v) -u d g A j(v,u) 

Fd,AA u > v ) = ~ d 2 • 17 

u a — v a 

This equation says that the 'crossing symmetry deficit' of all the fields in the OPE, normalized to 
the deficit of the unit operator, has to sum up to 1. 



Let us view Eq. (16) as a linear relation in the vector space of functions of two variables f(u, v). 
Then it can be given the following geometric interpretation. As we keep the CFT spectrum fixed 
and vary the squared OPE coefficients c\ l > 0, the vectors in the RHS fill in a convex cone 
generated by the functions F^^i- The sum rule says that the function f(u, v ) = 1 must belong to 
this cone (see Fig 3.4^ ). 



If we start imposing restrictions on the CFT spectrum, for example by demanding that there 
should be a gap in the scalar sector: A > A* for I = 0, this reduces the list of vectors generating 
the cone, and a fortiori the cone itself. It may well happen that the new reduced cone no longer 
contains the function / = 1, Fig |3.4b . A spectrum leading to such a cone cannot be realized in 
any CFT. 

If the situation in Fig |3.4b occurs, then, since the cone is convex, one can always find a 
hyperplane passing through the origin and separating / = 1 from the cone, Fig 3Alp . In analytical 



language, this means that there exists a linear functional A taking values of opposite sign on / = 1 
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Figure 3: Geometric interpretation of the sum rule: (a) the sum rule has a solution <^ / = 1 
belongs to the cone; (b) the assumed spectrum is such that the sum rule does not allow for a 
solution f = 1 does not belong to the cone; (c) in the latter situation, a hyperplane (the zero 
set of a linear functional) can be found separating f = 1 from the cone. 

and on the functions generating the cone: 

A[l] < 0, A[F dA>l ] > (18) 

In practice, the functional may be build up as a linear combination of the partial derivatives with 
respect to z and z at the democratic point z = z = 1/2. 




Figure 4: Geometric interpretation of Eq. (19). As t increases, the vector 1 — tF d ^ 2 eventually 
leaves the cone. 

So far we have described the method used in PQ,[2] to constrain the maximal allowed gap in 
the scalar sector. In order to constrain the size of the OPE coefficient c^, we proceed as follows 
[3j^J Let us rewrite the sum rule by transferring a part of the stress tensor contribution into the 
LHS: 

1 - tF dA2 = (c| 2 - t)F dA2 + £ c 2 A jF dA i (19) 

(A,0/(4,2) 



3 We choose 04,2 for definiteness; the method in fact works for any OPE coefficient. 
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The geometric interpretation of this equation is that the t-dependent vector l—tF d ^ 2 (u, v) belongs 
to the same cone as before as long as t < c\ 2 . In other words, the maximal allowed value of c\ 2 can 
be determined as the value t = t cr for which the curve 1—t F^p(u, v) crosses the cone boundary, 
Fig. |4j Analytically, we can detect that the crossing happened if there exists a linear functional 
such that 

A[F dA ,i] > (20) 

for all functions generating the cone, and 

A[l-tF dA2 ] = 0. (21) 

Note that in the present situation the function / = 1 must of course belong to the cone, otherwise 
the CFT simply does not exist and there is no point of discussing an upper bound on the OPE 



coefficients. Thus we are assuming from the start A[l] > 0, unlike in (18) 



Since the functional is linear, Eq. (21) is satisfied for 

* = A[l]/A[F d>4l 2], (22) 
and for larger t the functional will become negative as long as A [^,4,2] > 0. Thus we obtain the 



following result: each functional A satisfying (20) gives a bound on the maximal allowed value of 



r 2 ■ 
c 42 . 



max4 2 <A[l]/A[F d)4 , 2 ]. (23) 
This bound can be optimized by choosing the functional judiciously. 

The method just described was first applied in [3] to constrain the size of the OPE coefficients 



of scalar operators, while here we will use it to constrain the size of 04,2, which via ( 13 ) will give us 
a lower bound on Ct- Another difference from [3j is that we will study how the bound improves 
as a function of the assumed gap in the scalar sector of the OPE. 



4 Results 

We will now present our numerical results. First of all, let us consider the most general case 
when we are not making any assumption concerning the gap in the scalar sector of the OPE. This 
means that the scalar operators appearing in the OPE are allowed to have any dimension A > d. 
Operators with lower dimensions are a priori excluded if is the lowest dimension scalar. Under 
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this assumption, we use the method of linear functionals to bound c\ 2 from above. For this study, 
we choose linear functionals of the form 

A ^= £ ^W/ia=6=0 ( 24 ) 
n,m even, 
0<n+m<N 

z = 1/2 + a + b, z = l/2 + a-6. (25) 

As advertized, we are working around the democratic point z = z = 1/2. The fact that we are 
choosing z and z as real and independent can be interpreted as a Wick rotation to the Minkowski 
space pp. The functional only contains even derivatives because the functions F^a,i are even in 
both a and b pQ. 



We will choose Ao,o = 1 to have A[l] = 1. Then to optimize the bound (23), the coefficients of 
the functional must be chosen so that 

A[F dA2 ] ^max, (26) 



subject to the constraints (20), which in our case mean 



M F d,A,o\ > for all A >d, 

M F d,A,i) ^ for all A > / + 2, I = 2, 4, . . . (27) 

We will consider the functionals with the maximal derivative order up to N = 16.. Pushing 
to higher N values is likely to somewhat improve the bound. In principle N as large as 18 were 
demonstrated feasible in this kind of studies [2J. 



Eqs. (26), (27) define an optimization problem for the coefficients X m>n . The constraints are 



given by linear inequalities, and the cost function is also linear, which makes it a linear pro- 



gramming problem. Although the number of constraints in (27) is formally infinite, they can be 
reduced to a finite number by discretizing A and truncating at large A and I, where the con- 
straints approach a calculable asymptotic form pQ. The reduced problem can be efficiently solved 
by well-known numerical methods, such as the simplex method. A found solution can be then 
checked to see if it also solves the full problem. This procedure was developed and successfully 
used in PQ,0,[3]. 

Using this procedure, we computed a bound on c\ 2 from above, which via (13) translates into 



a bound on CV from below. The latter bound is plotted in Fig. [5] as a function of the dimension 
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Figure 5: The lower bound on the central charge Ct in terms of the dimension d of the lowest- 
dimension scalar primary. The stronger bound (upper blue curve) is obtained with N — 16. For 
comparison we give a weaker bound obtained with N = 12 (lower red curve), which corresponds 
to the horizontal axis A* = d in the following Fig. [#| The horizontal dashed line Ct = 4/3 shows 
where our bound stays above the free scalar central charge. 

of in the range 1 < d < 2. We plot our best bound for N = 16 and, for comparison, a weaker 
bound obtained with a smaller value N = 12. 

Postponing the discussion to the next Section, let us now consider what happens with the 
bound in presence of a gap in the scalar spectrum. In other words, we now assume that the 
first scalar operator in the 0x0 OPE has dimension A* strictly bigger than d. Technically, this 



problem is analyzed exactly as the previous one, except that the first set of constraints (27) is 
replaced by a shorter list: 

A[Fd, A ,o] > for all A > A* . (28) 

Because of considerable computer time involved, we solved this problem by using linear functionals 
with N = 12 only. The bound is given in Fig. [6] as a contour plot in the d, A* — d plane. On the 
horizontal axis A* = d the bound reduces to the N = 12 bound from Fig. |5j Naturally, when A* 
increases, the bound on Ct gets stronger. The white region in upper left corresponds to 

A, > 2 + 0.7(d - 1) 1/2 + 2.1(d - 1) + 0.43(d - 1) 3/2 (29) 

and is excluded, since such a large gap cannot be realized in any CFT according to the results of 
0. 
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A text file with the coefficients of linear functional used to derive the shown bounds is included 
together with this arxiv submission. 




d 

Figure 6: Contour plot of the Ct lower bound as a function of d and of the gap A* — d, where A* 
is the dimension of the first scalar in the 0x0 OPE. The gap is nonnegative, since we assume 
that is the lowest dimension scalar. On the horizontal axis the bound reduces to the N = 12 
curve in Fig. ^ The lighter green color marks the region where the bound is above C^ ee = 4/3, 
while in the darker red region the bound is below this value. As the gap increases, the bound gets 
stronger, so that a rather weak assumption about the gap is already enough to have Ct > Cfi ee . 

We end this section with a comment concerning the case of 2D conformal field theories. Recall 
that in PQ,[2] the maximal allowed gap in the scalar spectrum was studied for the 2D case in 
parallel with 4D. This was instructive since it allowed us to compare our bounds with the known 
OPEs in the 2D minimal models. The analysis is feasible because the 2D conformal blocks are 
known in a form just as simple as ^ (in odd dimensions similarly simple expressions are not 
available). Analogously, in the course of this project we have looked at the lower bounds on the 
Virasoro central charge c in the 2D CFTs, in the same d, A* plane (d > as appropriate in the 
2D case). We do not present them here because, in the range that we considered, the found 
lower bounds were smaller than 1. Since all unitary 2D CFTs with c < 1 are classified (these are 
precisely the unitary minimal models [17]). our bounds do not add any new information in this 
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case. 



5 Discussion 



Figs. 5]6 contain our advertized main results: universal lower bounds on the stress tensor central 



charge Ct- Fig. [5] gives a bound as a function of the dimension d of the lowest-dimension scalar 
present in the CFT. More precisely, the only requirement on is that the OPE 0x0 not contain 
any scalar of dimension less than d; this requirement is trivially satisfied if if lowest-dimension. 

The first interesting point about this bound is that in the limit d — > 1 it approaches the free 
scalar central charge value C^ ee = 4/3, see Eq. (J3J. In other words, our method shows that free 
theory limit is approached continuously. This is just as in previous work, where we proved that 
as d — > 1, the first scalar in the 0x0 OPE must have dimension below 2 [Tj,[2], and the 3-point 
function (000) must approach zero [3]. 

Next, we see that for 1 < d < 1.4 our bound stays above C^ ee , thus showing that an interacting 
theory necessarily has larger central charge than the free one. This is also rather interesting. 
Unfortunately, for larger d our bound drops below Cfi ee . We do not know if this means that there 
are CFTs with Ct < C^ ee . More likely, this indicates that our bound is not best-possible in this 
range. One could speculate that the best-possible bound should stay above C^ ee in the whole 
range 1 < d < 2. The fact that it should necessarily come down to C^ ee (or lower) for d = 2 can 
be inferred by considering the dimension 2 operator if 2 in the free scalar theory and its OPE with 
itself. 

Note that we could also derive a bound without using the assumption A > d, which would be 
applicable to any scalar, not just the lowest-dimension one. We have in fact derived also such a 
general bound, although we do not show it here. We found that this general bound differs little 
from the bound shown in Fig. [5] in the region of small d, say for d < 1.3. Thus this general bound 
could be useful if the lowest dimension scalar has dimension very close to 1, while the second- 
lowest is somewhat above 1. However, in the region of larger d, d > 1.7 -j- 1.8, the general bound 
drops to zero. This happens for the same technical reason that the bounds on the scalar OPE 
coefficients in [3] were blowing up around this value of the operator dimension. Because of this, 
in this paper we focused on the lowest dimension scalar, which allowed us to obtain a nontrivial 
bound in the full considered range of d. 

Now let us discuss Fig. [6j which gives the lower bound on Ct as a function of d and A*. 
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Here A* is the dimension of the lowest-dimension scalar in the OPE 0x0, assumed to be above 
d. Again, this assumption is trivially satisfied if is the lowest-dimension scalar present in the 



theory. Moreover, by the results of [T],[2] A* is limited from above by the bound given in Eq. (29). 

On the horizontal axis A* = d the bound in Fig. [6] reduces to the one shown in Fig. [5j while 
for larger A* it naturally gets stronger. In fact we see that A* somewhat bigger than d is already 
sufficient to raise the bound above C^ ee for all d (the lighter green region in the plot). The points 
with A* ~ 2d (i.e. with an approximate factorization of operator dimensions) belong to the green 
region by a big margin. 

In summary, we have shown in this work that if a unitary 4D CFT is non-trivial (in that it 
contains at least one primary scalar operator), then its central charge Cy cannot be arbitrarily 
low. We presented a universal bound on Ct as a function of the dimensions of the lowest and 
second-lowest scalar. We hope that these bounds will be helpful in future efforts to chart the 
"landscape" of 4D conformal theories. 

A relation like the one we derived, viewed from the AdS/CFT perspective (although of course 
we cannot do it since we are not at large N), would represent a lower bound on the Planck mass. 
One could then speculate that our result belongs to the same class of constraints on quantum 
field theory as the gravity as the weakest force conjecture [18]. Unfortunately our result cannot 
be directly applied to a phenomenon as fascinating and unavoidabe as gravity, but it has a non- 
negligible consolation that it follows from a rigorous mathematical analysis. 

We believe more general constraints of the type discussed in this paper lie ahead ready to be 
uncovered. 

Note added. The morning of the day this paper was submitted to arxiv, a nice paper 
appeared which, among other things, also derives lower bounds on the central charge. 
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